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数据 质量 控制 方法 与 效果 分 析 
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摘 要 : 针对 隆 德 站 地 基 微 波 辐 射 计 温度 资料 ,根据 数据 本 身 规律 及 各 类 疑 误 数据 表现 形式 ,结合 
历史 探 空 资料 ,利用 垂直 变化 强度 极 值 检查 、 标 准 差 检查 、 极 值 检 查 、 奇 异 值 检查 和 伪 值 检查 5 种 检 
查 方法 进行 质量 控制 实验 ,与 同期 探 空 资料 比 对 ,分 析 其 应 用 效果 。 结 果 表 明 : 各 质量 检查 结果 对 
数据 质量 均 有 一 定 的 区 分 能 力 , 且 对 控制 参数 是 敏感 的 ;数据 质量 以 晴空 最 优 、 云 天 次 之 、 降 水 入 
差 , 各 高 度 层 数据 质量 基本 相当 。 质 量 控制 结果 相关 性 分 析 表 明 各 高 度 层 上 晴空 相关 性 最 好 ,云天 
次 之 ,降水 稍 差 ,相关 性 都 处 在 较 高 的 水 平 上 ;降水 情况 下 质量 控制 效果 最 显著 。 
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随 着 气象 科学 事业 的 不 断 深入 与 发 展 ,新 型 仪 
器 仪表 有 力 地 推动 着 气象 观测 能 力 和 技术 手段 的 提 
高 ,气象 观测 的 数据 量 不 断 增加 ,然而 气象 资料 质量 
往往 受到 测 站 位 置 .观测 仪器 ,观测 技术 .观测 时 间 、 
观测 方法 等 非 气象 因素 的 影响 2 -31 。 气 象 资料 的 质 
量 控制 是 气象 资料 处 理 中 一 项 十 分 重要 的 工作 六 ， 
观测 研究 必须 对 其 所 用 的 资料 首先 进行 质量 检查 与 
处 理 , 这 样 才能 揭示 真实 的 天 气 .气候 变化 特征 与 规 
律 。 因 此 ,结合 气象 数据 特点 和 新 型 仪器 仪表 的 
言 息 获 取 原 理 与 方法 ,开展 新 型 气象 仪器 仪表 的 观 
测 数据 验证 和 质量 控制 研究 ,可 为 气象 观测 和 数据 
应 用 提供 技术 保障 9 。 

地 基 微 波 辐射 计 是 观测 不 同 高 度 大 气 参 数 的 新 
型 探测 设备 ,可 同时 反 演 高 垂直 分 辨 率 的 大 气温 度 、 
湿度 液态 水 廓 线 和 大 气 垂直 积分 水 汽 . 云 液态 水 总 
量 等 数据 ,具有 可 无 人 值守 连续 工作 ,高 时 间 和 空间 
分 辨 率 .操作 简洁 方便 等 优点 ,是 对 常规 探 空 的 有 益 
补充 ,已 逐渐 成 为 遥感 大 气温 度 ,湿度 .液态 水 廓 线 
以 及 大 气 水 汽 和 云 液态 水 总 量 的 有 力 工 具 " "|。 
利用 地 基 微 波 辐射 计 , 可 高 效 .准确 观测 降水 前 后 云 
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中 水 汽 液态 水 含量 分 布 和 变化 ,为 暴雨 天 气 预 报 、 
AT HAR (EEE HEAT Fa WE BPP RAE 
节 和 天 气 条 件 下 微波 辐射 计 探 测 资料 的 精度 有 差 
FeO) ,微波 辐射 计 反 演 系数 不 当 、 天 线 日 上 的 积 
水 .电磁 干扰 .电源 .通信 等 外 界 影响 会 导致 观测 数 
据 出 现 异 常 , 因 此 ,开展 系统 的 质量 控制 是 合理 有 效 
使 用 微波 辐射 计 探测 资料 的 基础 。 

针对 常规 气象 观测 数据 ,国内 外 气象 工作 者 提 
出 了 多 种 质量 控制 方案 。COLLINS: 给 出 了 针对 
大 气 探 空 资料 中 高 度 和 温度 的 质量 控制 方法 ,发 展 
和 完善 了 常规 质量 控制 方法 。DURRE'” 针对 全 球 
探 空 数 据 集 设计 了 一 套 完整 的 质量 控制 系统 ,对 温 
度 超过 气候 极 值 垂直 一 致 性 差 . 伪 值 .时 间 突 变 等 
多 种 类 型 的 疑 误 数据 设计 了 不 同 的 检查 方法 。 能 安 
元 ' 沾 介绍 了 北欧 国家 对 实时 和 非 实时 气象 资料 进 
行 质量 控制 的 流程 以 及 采用 的 方法 .技术 等 。 张 志 
富 ' 研 制 了 适用 于 自动 站 土壤 水 分 小 时 数据 的 质 
量 控制 方案 。 

微波 辐射 产品 属 非常 规 气 象 探测 数据 ,产品 一 
般 分 为 3 级 ,其 中 Level-0 为 原始 电压 ,Level-1 为 亮 
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温 ,Level-2 为 温 .湿度 廓 线 及 垂直 积分 水 汽 .液态 水 
含量 等 产品 。 因 其 时 间 分 辩 率 更 高 ,探测 原理 .错误 
来 源 . 疑 误 数据 表现 形式 与 常规 气象 资料 存在 差异 ， 
针对 常规 气象 观测 资料 的 质量 控制 方法 不 完全 适用 
于 微波 辐射 计 习 ; 。 鉴 于 此 ,国内 部 分 气象 工作 者 对 
微波 辐射 计数 据 质量 控制 进行 了 研究 。 李 青 '“ 等 
基于 辐射 传输 理论 ,检验 了 一 台 多 通道 微波 辐射 计 
亮 温 数据 的 “晴空 ”观测 和 数值 模拟 结果 之 间 的 一 
致 性 。 朱 雅 航 等 2 针对 地 基 微 波 辐射 计 , 研 究 其 亮 
温 数据 的 质量 控制 方案 ,提出 极 值 . 时 间 一 致 性 、 辆 
射 传输 计算 .多 通道 亮 温 交 又 检查 。 歼 雪 等 "通过 
个 例 分 析 微 波 辐射 计 观 测 值 与 计算 值 的 一 致 性 ,并 
判断 各 通道 观测 值 的 合理 性 。 这 些 质量 控制 方法 主 
要 针对 Level-1 数据 展开 ,而 针对 Level-2 数据 研究 
较 少 。 目前 ,国内 应 用 较 多 的 是 以 美国 MP-3000 A 
为 代表 ,以 " 串 行 常规 合成 器 变频 技术 ”为 基础 的 地 
基 微 波 辐射 计 。 而 德国 RPG 公司 研制 的 RPG- 
HATPRO-G4 型 地 基 多 通道 微波 辐射 计 ( 以 下 简称 
RPG ) 采 用 “多 通道 并 行 测 量 技 术 ”, 其 工作 原理 与 
MP-3000 A 型 地 基 微 波 辐射 计 不 同 , 有 其 自身 的 技 
术 特 点 i。 本 研究 利用 甘肃 平凉 气象 站 历史 探 空 
资料 和 宁夏 隆 德 县 气象 站 RPG-HATPRO-G4 型 地 基 
多 通道 微波 辐射 计 资料 ,在 总 结 前 人 数据 质量 控制 
方法 的 基础 上 ,根据 数据 本 号 的 特点 及 各 种 类 型 疑 
误 数据 表现 形式 ,首次 给 出 了 一 套 针 对 RPG- 
HATPRO-G4 型 地 基 多 通道 微波 辐射 计 Level-2 温度 
数据 的 质量 控制 方法 ,为 业务 中 有 效 使 用 其 温度 资 
料 提供 技术 支持 。 


1 数据 与 方法 


数据 处 理 
选用 隆 德 气象 站 (35. 37° N, 106. O7°E, He HK 
2 078 m)2017 年 6 月 10 日 ~2017 年 12 月 31 日 
(仪器 于 6 月 9 日 完成 标定 ) RPG 温度 数据 展开 研 
究 , 其 时 间 分 辩 率 为 1 min ,垂直 高 度 (五 )0 ~ 10 km 
内 共 93 层 ,垂直 分 辩 率 见 表 1 。 

选用 甘肃 平凉 站 (35.55 °“N,106. 67 “上 ,海拔 
1 468 m, 有 中 隆 德 气 象 站 49 km)2013 一 2017 年 6 ~ 12 
月 每 日 两 次 (08 时 ,20 时 ) 的 探 空 数据 ,作为 确定 质 
量 控制 净值 的 样本 。 每 组 数据 时 间 分 辩 率 为 1 s, E 
直 分 状 率 为 6 ~7 m。 探 空 资料 已 通过 台 站 级 质量 
控制 ,可 直接 使 用 。 


1.1 


R1 RPG 温度 廓 线 垂直 分 辨 率 / m 
Tab.1 Vertical resolution of the RPG temperature 


profile / m 
HO è ® ä 率 a 率 
0~10 10 500 ~ 1 200 40 3 500 ~4 500 160 
10 ~25 15 1 200 ~ 1 800 60 4500 ~6 000 200 
25 ~ 100 25 1 800 ~2 500 90 6 000 ~9 800 300 
100 ~500 30 2 500 ~3 500 120 | 9 800 ~10 000 200 


对 RPG 数据 和 探 空 数据 做 如 下 处 理 : 

(1) 统一 高 度 :将 探 空 数据 插值 到 对 应 RPG93 
个 高 度 层 上 , 共 获 得 探 空 廓 线 样本 2 139 个 , 探 空 数 
据 样 本 197 065 个 。 

(2) 统一 时 间 : 将 RPG 时 间 格 式 转换 为 北京 
时 ,根据 2017 年 6 H ~ 2017 年 12 月 插值 后 探 空 数 
据 起 止 时 间 ,选取 对 应 时 间 范 围 内 的 RPG 温度 数 
HE , 共 获 得 温度 廓 线 样本 9 679 个 ,数据 样本 900 147 


个 
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1.2 质量 控制 方法 
针对 RPG 温度 廓 线 数 据 特 点 ,根据 ZAHUMEN- 
SKY") hyper) 、 王 伯 民 于 等 质量 控制 设计 思路 和 
方法 ,设计 垂直 一 致 性 检查 . 极 值 检查 和 时 间 一 致 性 
检查 3 类 检查 方法 ( 含 5 个 单项 检查 方法 ) ,并 最 终 
给 出 组 合 和 加 权 2 种 综合 检查 结果 (图 1)。 综 合 考 
虑 检查 结果 的 区 分 度 和 效果 , 均 将 3 类 .5 项 检查 结 
果 分 成 4 档 , 分 别 用 数据 质量 标识 符 0.3 6 9 区 分 
标识 1。 
1.2.1 垂直 一 致 性 检查 ” 相 邻 高 度 上 温度 差异 应 
在 一 定 范围 内 , 探 空 温 度 资料 质量 控制 一 般 都 包含 
垂直 一 致 性 检查 '”) 。 常 见 不 满 足 垂直 一 致 性 的 情 
况 主要 有 年 直方 向 存在 奇异 值 . 层 际 变化 大 等 特 
fe! =, RPG 温度 廓 线 数据 记 为 7, 垂直 变化 强度 记 
J dT, (HA:C + (100 m) 7!) ,计算 公式 如 下 : 
-T, 


(1) 


式 中 ;7 Tra A h h +1 AWE; H, Hp, AT 
的 高 度 。 
(1) 垂直 变化 强度 极 值 检 查 

垂直 变化 强度 极 值 检 查 就 是 dT, 应 在 历史 探 空 
样本 垂直 变化 强度 (ds ) 的 最 大 、 最 小 值 之 间 取 值 。 
考虑 dT, 存在 显著 的 季节 性 , 且 在 各 高 度 上 存在 差 
异 ,开展 检查 时 ,以 自然 月 各 高 度 层 为 划分 单位 , 参 
考 探 空 资料 统计 结果 ,针对 6 ~12 月 92 层 分 别 设置 


tó Wwy 


垂直 变化 强 垂直 变化 强度 
度 极 值 检查 标准 差 检查 
(n,) (n,) 


垂直 一 致 性 综合 
检查 (N ) 


探 空 温度 (Ts) 


RPG 温 度 (了 T) 


设置 质量 控制 符 门 限 ) 


质量 控制 QC HERE, 


奇异 值 检查 EERE 
(n,) (n,) 


时 间 一 致 性 综合 
检查 (N,) 


图 1 RPG 温度 数据 质量 控制 方案 流程 


Fig.1 Flow chart for quality control procedure of RPG temperature data 


检查 冰 值 ,最 大 最 小 值 分 别 记 为 ATS pay SETS pin 9 FE 
义 6 为 二 者 之 差 除 以 100。 依 据 所 获 阔 值 , 逐 廓 线 
逐 层 开展 检查 ,数据 质量 标识 符 用 n, 表示 ,含义 见 
表 2。 
(2) 垂直 变化 强度 标准 差 检查 

垂直 变化 强度 标准 差 检查 以 温度 层 际 变化 不 应 
太 大 为 依据 ,将 垂直 变化 强度 标准 差 记 为 o, 
(2) 式 ,其 中 m 为 廊 线 样本 数 。 


o= JEL (dn) 0) 


垂直 变化 强度 标准 差 检 查 采 用 动态 国 值 , 即 当 
日 当 层 o 值 由 当天 为 中 心 , 步 长 为 15 d 的 该 层 dy 
样本 求 得 。 根 据 朱 雅 笠 分 档 和 质量 标识 方法 , 逐 
廓 线 逐 层 进行 检查 ,对 满足 一 定 条 件 的 数据 进行 质 
量 标识 。 考 虑 六 盘山 区 相 邻 时 刻 温 度 变 化 及 垂直 变 
化 强度 较 东 部 治 海地 区 大 ,质量 控制 国 值 进行 一 定 
尺度 放宽 ,数据 质量 标识 符 用 n, 表示 ,含义 见 表 3。 


表 2 垂直 变化 强度 极 值 检查 数据 质量 标识 符 n, 取 值 解释 
Tab.2 Values of quality identifier n, for extreme value 
check of vertical change intensity and the 


corresponding meaning 


ny 检查 条 件 

0 dT, €e [dTs ymin ydTS mar | 

3 dT, € [dTs nmin 38, dTs hmin | U [dTs pmax 1478 hmax +38] 

6 dT, € LATS jpn 68, dTs jin — 35] U [dTs may +38, dTs pmax +68] 
9 dT, e[ -œ ,dTs,,m 68] U [dTs mas +65, + | 


表 3 垂直 变化 强度 标准 差 检查 数据 质 
量 标识 符 n, 取 值 解释 
Tab.3 Values of quality identifier n, for standard 
deviation check of vertical change intensity and the 


corresponding meaning 


Ny 检查 条 件 
0 dT, e [dT, +20] 


3 dT, e [dT, -2.50 ,dT, -20] U [dT, +20,dT, +2.50] 


6  dT,eldT,-30,dT,-2.50]U[dT, +2.50,dT, +30] 


\o 


dT,e[ -œ ,dT, -30] U[dT, +30, +% ] 


(3) 垂直 一 致 性 综合 检查 

结合 加 权 判 别 技术 和 数据 质量 分 档 标 识 方 
YEON ,垂直 一 致 性 综合 检查 数据 质量 标识 符 用 N 
表示 ,含义 见 表 4。 
1.2.2 meee 极 值 检查 就 是 7 应 在 历史 探 空 
样本 温度 最 大 、 最 小 值 之 间 取 值 。 开 展 极 值 检查 时 ， 
以 自然 月 各 高 度 层 为 划分 单位 ,参考 探 空 历史 资料 
统计 结果 ,针对 6 月 ~12 月 93 层 分 别 设置 极 值 检 
AE BOE PDI Ts na FU TS pin ,定义 i 为 二 者 之 差 除 
以 100 ,数据 质量 标识 用 N, 表示 ,含义 见 表 5。 


表 4 垂直 一 致 性 综合 检查 数据 质量 标识 符 N, 取 值 解释 


Tab.4 Values of quality identifier N, for comprehensive 


check of vertical consistency and the corresponding meaning 


Ni 检查 条 件 N, 检查 条 件 
0 n +m =0 6 ny +m =6 
3 ny +m =3 9 ny +m 二 9 
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R5 极 值 检查 数据 质量 标识 符 N, 取 值 解释 
Tab.5 Values of quality identifier N, for 


extreme value check and the corresponding meaning 


Ny 检查 条 件 

0 Te [TS min» TSmax] 

3 Te [Ts min 3t, Ts min | U [ Ts max s Smar +3t] 

6 T ETS pin 6t, Ismin 3t] OL Ts mar +3t, TS max +6t] 
9 Tel -œ ,Tsmn 6t] U[Tsma +6t, + | 


1.2.3 时 间 一 致 性 检查 

(1) 奇异 值 检查 

奇异 值 检 查 ,适用 于 温度 随时 间 变 化 存在 突变 
的 数据 。 定 义 相 邻 时 次 变化 值 超过 其 平均 值 m 倍 
标准 差 的 数据 为 疑 误 数据 (定义 见 (3) (4) 式 ) 。 


a ah, 


>mxo (3) 


dT, =T, -T,; (4) 


式 中 :d7, 为 温度 随时 间 的 变化 值 ;7, 和 也 ,分别 为 
it 和 (t+1) 时 刻 温度 ;m 为 偏离 标准 差 的 倍数 ;o 定 
义 见 (2) 式 。 质 量 控制 国 值 o 仍 采用 动态 阔 值 , 数 
据 质量 标识 符 用 ny 表示 ,含义 见 表 6。 

(2) 僵 值 检查 

僵 值 检查 适用 于 温度 随时 间 变 化 而 无 变化 的 数 
据 , 依 据 相 邻 时 次 廊 线 值 连续 不 变 的 条 数 ( 记 为 m) 划 
分 数据 质量 等 级 ,质量 标识 符 用 ns 表示 ,含义 见 表 7。 

(3) 时 间 一 致 性 综合 检查 

时 间 一 致 性 综合 检查 数据 质量 标识 符 用 M, 表 
示 , 含 义 见 表 8。 


表 6 奇异 值 检查 数据 质量 标识 符 n 取 值 解释 
Tab.6 Values of quality identifier n, for singular 
value check and the corresponding meaning 
na 检查 条 件 
dT, e [dT, +20] 

dT, e [dT, -30 ,dT, -20 ] U[ dT, +20,dT, +30] 
dT, e [dT, -40 ,dT, -30 ] U [ dT, +30 ,dT, +40] 

dT,e[ -œ ,dT, -4¢]U[dT, +40, +% ] 


O o uU oO 


R7 僵 值 检查 数据 质量 标识 符 ns 取 值 解释 
Tab.7 Values of quality identifier n, for invariable 


value check and the corresponding meaning 


0 m=0 6 me [6,9] 
3 me [2,5] 9 me[10, +æ] 


表 8 时 间 一 致 性 综合 检查 数据 质量 标识 符 N, 取 值 解释 
Tab.8 Values of quality identifier N, for comprehensive 


check of time consistency and the corresponding meaning 


Ny 检查 条 件 N; 检查 条 件 
0 ng +n; =0 6 ng +n; =6 
3 n4, +n; =3 9 n, +n; >=9 


1.2.4 综合 检查 “完成 3 类 检查 后 ,进行 综合 检 
查 。 综 合 检 查 中 ,组 合 检 查 定义 为 3 类 检查 质量 标 
识 符 的 各 类 组 合 (N = N NmN: ) , N 各 数位 值 越 大 , 相 
应 数据 疑点 越 大 ;加 权 检 查 定义 为 3 类 检查 质量 标 
识 符 之 和 (下 =N +N, + NN,) ,下 值 越 大 ,数据 疑点 越 
大 。 


2 结果 与 分 析 


2.1 三 类 质量 检查 结果 分 析 

为 对 比 不 同 天 气 背景 ` 不 同 高 度 的 数据 质量 检 
查 结果 ,将 N,N, AN, 的 93 层 数 据 间隔 500 m 分 
FA 20 层 , 并 按 不 同 天 气 背 景 条 件 分 别 进 行 频 率 分 
布 统计 (图 2)。 其 中 晴空 (各 高 度 层 相对 湿度 < 
85%) 云天 (无 降水 且 有 1 层 及 以 上 相对 湿度 
=85% ) 、 降 水 样本 数 分 别 为 477 834、326 802、 
95 511 个 。 晴 空 背 景 下 各 高 度 层 N, = 0、N, = 0、 
N=0 的 数据 占 比 均 分 别 高 于 94. 28% 、99. 91% 、 
92. 64% ;云天 背景 下 各 高 度 层 N, =0.N, =0 Ns =0 
的 数据 占 比 均 分 别 高 于 95. 53%、98. 26% 、 
93.27% ;降水 背景 下 各 高 度 层 N, =0 WN, =0.N, =0 
的 数据 占 比 均 分 别 高 于 91. 83% 、91. 62% 、 
78.77% 。3 类 检查 方法 对 数据 质量 均 有 一 定 的 指 
示 能 力 ,数据 质量 均 以 晴空 最 优 、 云 天 次 之 .降水 最 
低 ,N =9 N, =9 N; =9 均 以 降水 背景 贡献 最 大 。 
2.2 综合 质量 检查 结果 与 区 分 度 分 析 

综合 检查 中 各 位 数 的 取 值 或 下 值 的 大 小 ,对 
应 着 一 定 的 数据 质量 区 分 能 力 , 称 为 区 分 度 '”  。 对 
N 进行 全 样本 统计 ,对 下 进行 不 同 天 气 背 景 及 全 样 
本 的 分 类 统计 ,具体 分 析 如 下 。 

(1) 全 样本 900 147 个 数据 随 N 的 变化 如 表 9 
所 示 ,N=000 数据 个 数 为 811 859, 占 数据 总 量 的 
90. 19% ;NN 各 位 数 最 大 的 组 合 为 V =399 ,其 数据 个 
数 仅 为 94, 占 数据 总 量 的 0.01% 。 

(2) 全 样本 及 不 同 天 气 背 景 下 数据 样本 随 F 
的 变化 如 表 10 所 示 , 取 <3, 则 全 样本 、 晴 空 、 云 
天 、 降 水 情况 下 ,可 能 高 质量 数据 分 别 占 各 自 数 据 量 
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x 2 
tid wN 
Mo MSs Msc Do 
a) 晴空 N (b,) ZN, (c) 晴空 N, 
100 100 100 
95 
90 
95 95 
85 
80 
90 90 75 
(a)BRN, (b,) BEN, (C) BEN, 
100 100 100 
x 95 
w 90 
g 9 95 
= 85 
BE 80 
90 90 75 
(a,) 降水 NN (b,) EKN, (cs) JKN, 
100 100 100 
95 
90 
95 95 
85 
80 
90 90 75 


0 1 2 3 4 5 6 7 8 39) 10 


OL2s 4 5 6 7 8 9 10 
BE /km 


OL23 4 5 6 7 & 9 10 


图 2 不 同 天 气 背景 ,不 同 高 度 ,不 同 检查 方式 微波 辐射 计 温 度数 据 质量 检查 结果 频率 分 布 


Fig.2 Frequency distribution of microwave radiometer temperature data quality check results in different weather conditions , 


at different heights and through different check methods 


RI 全 样本 组 合 检查 和 (NiN,N;) 对 应 数据 个 数 
Tab.9 Corresponding data number of combination 
check result N (N,N,N, ) with full sample size 
N 数据 个 数 N 数据 个 数 N 数据 个 数 N 数据 个 数 
000 811859 300 21592 600 6169 900 1 333 

003 34846 303 1197 603 295 903 61 


006 7390 306 333 606 187 906 15 
009 11729 309 475 609 144 909 4 
030 408 330 2 630 0 930 1 
033 66 333 0 633 0 933 0 
036 27 336 7 636 0 936 0 
039 121 339 13 639 0 939 0 
060 106 360 4 660 0 960 0 
063 36 363 1 663 0 963 0 
066 33 366 6 666 0 966 0 
069 135 369 13 669 0 969 0 
090 319 390 22 690 2 990 1 
093 327 393 16 693 1 993 0 
096 78 396 2 696 0 996 0 
099 677 399 94 699 0 999 0 


的 96.51% .97. 86% .97. 83% 和 84. 21% 。 数 据 个 
数 不 为 0 的 下 最 大 值 为 21 ,全 部 出 现在 降水 数据 
中 。 可 见 ,晴空 和 云天 背景 数据 质量 都 很 高 , 均 优 于 
降水 背景 数据 ,但 即使 是 数据 质量 相对 较 差 的 降水 


背景 数据 ,F<3 数据 个 数 依 然 占 降水 背景 数据 样本 
总 量 的 84.21% ,这 说 明 微波 辐射 计 观 测 性 能 优 恨 。 
2.3 敏感 度 分 析 

垂直 变化 强度 标准 差 检 查 结果 n, 和 时 间 一 致 
性 奇异 值 检 查 结果 n 依赖 于 各 自 的 e 值 。 如 果 分 
别 以 各 自 o 为 阅 值 参考 值 ,分 别 做 “放宽 闷 值 ”( 增 
K o) , 则 相应 地 “ns HUME” “n, 取 小 值 "的 发 生 
频率 应 该 增加 。 这 就 是 质量 检查 结果 对 质量 检查 阐 
BURE. n, 和 nm 分 别 以 各 自 o 作为 参考 
值 ,将 其 百分比 增 量 记 为 A,A 从 -25% ~ 25% 每 隔 
2.5% 取 一 值 ,重新 对 微波 辐射 计 温度 全 样本 数据 进 
行 检查 结果 的 频率 分 布 统计 ,如 图 3 所 示 。 

可 见 , 随 A 增 大 ,n, =0 和 n=0 的 频率 逐渐 增 
大 ,mm >0 和 ns >0 各 档 的 频率 均 逐 渐 减 小 ,表明 数 
据 质量 控制 结果 对 控制 参数 是 敏感 的 。 在 微波 辐射 
计 温 度数 据 的 后 续 使 用 中 ,用 户 可 根据 需要 ,通过 设 
置 质量 控制 符 门限 来 控制 使 用 的 数据 。 
2.4 相关 性 分 析 

根据 同期 探 空 资 料 起 止 时 间 选 取 获 得 RPG 温 
度 廊 线 样本 9 679 条 ,其 中 晴空 云天、 降水 样本 分 
别 为 5 138 、3514、1027 条 。 基 于 综合 质量 加 权 检 
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表 10 不 同 天 气 背 景 下 加 权 检 查 (Ni +N, +N, ) 对 应 数据 个 数 及 频率 
Tab.10 Corresponding data number and frequency of weighted check result F (N, +N, +N,) 


in different weather conditions 


p 全 样本 晴空 云天 降水 
数据 个 数 (频率 / % ) 数据 个 数 (频率 / % ) 数据 个 数 ( 频 率 / % ) 数据 个 数 (频率 / % ) 
0 811 859(90. 192) 436 926(91.439) 301 201(92. 166) 73 732(77. 197) 
3 56 846(6. 315) 30 681(6.424) 18 512(5. 665) 7 653(8.013) 
6 14 930(1.659) 7 255(1.518) 5 139(1.573) 2 536(2. 655 ) 
9 14 076(1. 564) 2 791(0. 584) 1 855(0. 568) 9 430(9. 873 ) 
12 1 235(0. 137) 161(0. 034) 83(0.025) 991(1.038) 
15 409(0.045) 17(0. 004) 12(0.004) 380(0.398) 
18 698(0.078) 3(0. 001) 695 (0. 728 ) 
21 94(0. 010) 94 (0. 098 ) 
100 p 10 
(a) n,=0, n,=0; (b)n,>0, n,>0 
98 上 
96 
x 94 s 
At vl 
R 92 长 
90 
88 | 
86 上 
85 1 上 L L L J 
-25 -15 5 0 5 15 25 25 
A/% 
图 3 垂直 变化 强度 标准 差 检查 (m% ) 和 时 间 一 致 性 奇异 值 检查 (mm ) 各 档 取 值 频率 随 A 的 变化 
Fig.3 Variation of frequency distribution at all data quality levels along with the change of A values in the cases of standard 


deviation check (n, ) of vertical change intensity and time consistency check (n,) of singular value check 


一 一 质量 控制 后 --- 质量 控制 前 


i (a) 晴空 107 NN. (ob) 云天 
9F 9 
8r 8 
Tr HA 
g 8) BEO £ 
z 5 上 z 5 z 
4 F 4 
3 上 3 
2 上 2 y 
lr 1 
0 : ， ， 0 : : i 
0.7 0.8 0.9 1.0 0.7 0.8 0.9 1.0 
相关 系数 相关 系数 
图 4 不 同 天 气 背景 下 加 权 检 查 ( 忆 =0 作为 质量 控制 符 门限 ) 质量 探 人 


© 
1 


N cp A CI for) | Oo No} 
T 


0.8 


相关 系数 


0.9 1.0 


ill Fd Ja AG AA BCR) 随 高 度 的 变化 


Fig.4 Variation of correlation coefficients (R) along with height before and after the quality control of weighted 


check method ( F) in different weather conditions 


tó wy 


Aiii F KEMENA TRAN, ECARD AS E] 
BY F (EIX RPG 温度 廓 线 数据 进行 质量 控制 ,再 
将 单 次 探 空 时 段 内 RPG 温度 资料 逐 层 进行 平均 , 共 
获得 与 探 空 温度 廊 线 一 一 对 应 的 廊 线 分 别 为 232、 
168 .64 条。 最 后 按 不 同 高 度 与 对 应 的 探 空 数据 进 
行 相关 统计 ,得 到 不 同 天 气 背 景 下 加 权 检 查 质 量 控 
制 前 后 相关 系数 随 高 度 的 变化 结果 ,如 图 4 所 示 , 图 
中 各 相关 系数 均 通 过 了 显著 性 水 平 为 0.01 的 显著 
性 检验 。 

可 见 , 各 高 度 层 上 晴空 背景 下 相关 性 最 好 ,云天 
次 之 ,降水 较 晴 空 和 云天 稍 差 , 但 总 体 相关 性 处 在 较 
高 水 平 上 。 以 下 =0 作为 质量 控制 符 门 限 为 例 , 质 
量 控制 后 对 应 整 层 相关 系数 平均 值 分 别 为 晴空 
0. 967 .云天 0. 934 降水 0.880。 降 水 情况 下 质量 控 
制 后 相关 系数 有 一 定 提高 ,=0、F S 3 所 对 应 数据 
的 相关 系数 与 质量 控制 前 数据 的 相关 系数 相 比 分 别 
提高 0.015 和 0. 008。 


3 结论 


综合 以 上 分 析 , 得 出 主要 结论 如 下 : 

(1) 各 类 数据 质量 检查 结果 对 数据 质量 均 具 有 
一 定 的 区 分 能 力 ,数据 质量 以 晴空 最 优 、 云 天 次 之 、 
降水 稍 差 ,微波 辐射 计 观 测 性 能 优良 。 

(2) 数据 质量 控制 结果 对 控制 参数 是 敏感 的 ， 
在 微波 辐射 计 温 度数 据 的 后 续 使 用 中 ,通过 设置 质 
量 控 制 符 门限 来 控制 使 用 的 数据 。 

(3) 质量 控制 前 后 ,各 高 度 层 上 晴空 背景 下 数 
据 相 关 性 最 好 ,云天 次 之 ,降水 较 晴 空 和 云天 稍 差 ， 
相关 性 都 处 在 较 高 的 水 平 上 ,降水 情况 下 的 数据 质 
量 控 制 效果 最 显著 。 
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Abstract: The quality control of meteorological data is a very important work in the processing of meteorological 
data , because the microwave radiation products are unconventional meteorological observation data ,so the method of 
quality control for conventional meteorological observation data is not fully applicable to microwave radiometer. In 
order to study the quality control method on the temperature data by RPG-HATPRO-G4 type ground-based multi- 
channel microwave radiometer and provide technical support for effective use of temperature data, this paper propo- 
ses five methods according to the rules of data itself and various types of the wrong data using the temperature data 
collected by RPG-HATPRO-G4 type ground-based multichannel microwave radiometer at Longde Station , Ningxia 
Province ,China, combined with the 5-year historical sounding data from Pingliang Station, Gansu Province , China 
which is 49 kilometers far from the Longde Station. The five methods include extreme value check of vertical change 
intensity , standard deviation check of vertical change intensity , extreme value check , singular value check and invar- 
iable value check. They are further combined into three types: comprehensive check of vertical consistency , extreme 
value check and comprehensive check of time consistency. And then the combination check and weighted check re- 
sults are given by synthesizing the three types of check results. Finally ,in order to analyze the effect of these meth- 
ods ,we applied these methods to the collected data and compared the results with the sounding data in the same pe- 
riod. The results show that the data quality in the case of clear weather is better than that in the case of cloud and 
precipitation weather, and the data in the cloud weather is better than that in the precipitation weather. For the 
abovementioned 3 weather conditions ,the data portion of high quality accounts for 97.86% ,97.83% and 84.21% 
of their corresponding overall data respectively. The data quality at each layer is basically equal, the data of poor 
quality are mostly contributed by the precipitation weather. Combination check and weighted check both have better 
ability to measure data quality and the respective quality control results are sensitive to the corresponding control 
parameters. The correlation analysis of the quality control results shows that the correlation is the best in the case of 
clear weather at each altitude ,followed by the cloudy weather,and the worst for the precipitation weather , and their 
correlation coefficients are all at high levels. The effect of data quality control is the most obvious in the case of pre- 
cipitation weather. The experiment in this study demonstrated this scheme of data quality control is feasible and pro- 
vided the mechanism to distill the temperature data from the RPG-HATPRO-G4 type ground-based multichannel mi- 
crowave radiometer for other applications. 


Key words: ground-based microwave radiometer ; temperature ; sounding data; quality control 


